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Abstract. We provide a new analytical approach to operator splitting for 
equations of the type ut = Au + B{u) where A is a linear operator and B 
is quadratic. A particular example is the Korteweg-de Vries (KdV) equation 
ut — uux +Uxxx = 0. We show that the Godunov and Strang splitting methods 
converge with the expected rates if the initial data are sufficiently regular. 

o 
o 

(N 

^ 1. Introduction 

The ubiquitous Korteweg-de Vries (KdV) equation 

^ Ut- UU^ + Uxxx = 

offers the perfect blend of the simplest nonlinear convective term uux and the 

p ^ simplest dispersive term u^^x'- The well-known smooth soliton solutions of the KdV 

equation interact in an almost linear fashion apart from a phase shift. Furthermore, 
. the KdV equation is completely integrable with an infinite family of conserved 

i^H quantities. This is a result of a subtle interaction between the Burgers term uUx 

and the Airy term Uxxx, as it is a well-known fact that the nonlinear Burgers 
^ equation Ut — uu^ — generically develops shocks in finite time while the linear 

I— I Airy equation ut + Uxxx — preserves all Sobolev norms. 

The initial value problem for the KdV equation with u|t=o = wo G -ff^, either 
^ on the whole real line or in the periodic case, has been extensively studied. An 

incomplete list of references is [2l|5]. For further discussions we refer to [3 [6]. 
The method of operator splitting, also called the fractional steps method, remains 
0^ a very popular method both for analysis and numerical computations of partial 

differential equations. Instead of including a long list of references to relevant works, 
\G we entrust the reader instead with [3] and the overview of the field given therein. 

However, we refer to [7] for rigorous analysis of the splitting method applied to the 
Schrodinger-Poisson and the cubic nonlinear Schrodinger equations. 

Formally operator splitting can be explained as follows. Let u{t) — <i>c(t)uo = 
<I>c(t;uo) € X, where X is some normed space, denote the solution of a given 
differential equation 

ut^C{u), te[0,T], u|t=o = ""o- 

Here C will typically be a differential operator in the spatial variable. Assuming 
that we can write C — A + B in a, natural way, the idea of operator splitting is that 

(1) u{t^) « ($B(At)$^(At))"wo, t„ = nAt, 
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where At ^ 1 . In the case of hnear ordinary differential equations this goes back 
to Sophus Lie. One of the reasons for its popularity is that operator splitting allows 
for a separate treatment of the equations Ut = A(u) and Ut = B{u); in particular 
this applies to the use of dedicated special numerical techniques for each of the 
equations, again we refer to [1] for a long list of examples and relevant references. 

In the context of the KdV equation, the first use of operator splitting was re- 
ported in the brief paper by Tappert ^ where it was applied as a numerical method. 
Apparently the first rigorous results appeared in [3] where a Lax-WendrofF result 
was proved: If operator splitting converges to some limit function, then the limit 
function is a weak solution of the KdV equation. In addition a systematic study 
of operator splitting as a numerical method was undertaken for the KdV equa- 
tion. More extensive rigorous results, specifically the convergence of the splitting 
approximations in a suitable functional space, were hampered by the apparent in- 
compatibility between the Burgers equation and the Airy equation. 

In the present paper we offer a new analytical approach to operator splitting for 
the KdV equation that will lead to rigorous convergence results (error estimates). 
Compared to earlier attempts, two new ingredients enter the present approach. 
First of all we actively use that the solution of the KdV equation remains bounded 
in a Sobolev space, that is, if uq € then remains bounded for 

t S [0,T]. This together with a bootstrap argument is used to secure the existence 
of a uniform choice of time step At that prevents the solution from any "Burgers" 
step from blowing up. Indeed the main problem in this approach is that the Airy 
equation produces small oscillatory waves that, when used as initial data for the 
Burgers equation, produce shocks. Secondly, since the splitting approximations are 
merely defined at the discrete times i = t„, to facilitate the convergence analysis we 
introduce an extension which is defined for all t G [0, T] . Concretely, we introduce an 
extension v which depends on an additional time variable t S [0, T], i.e., v — v(t, r), 
and let the evolution corresponding to each time variable be governed by one of the 
split operators in such a way that at each time level t — tn the extension w(t„, t„) 
coincides with the regular splitting approximation. This extension approach is 
different from the conventional one, where one lets "time run twice as fast" in 
each of the subintervals [t„, t„ -I- At/2] and [t„ -I- At/2, t„+i] (cf. the discussion and 
references in W). 

Formally the operator splitting ([T|) , called sequential or Godunov splitting, yields 
a first order approximation in At, that is, 

||u(t„) - ($A(At)$B(At))"uo|| = 0(At), t„ -> t, At -> 0, 
in an appropriate norm. We show that this holds rigorously for the KdV equation. 



More precisely, we prove in Theorem 2.4 that for mo G with s > 5 we have 

that for At sufficiently small 

||z;(t,t)-u(t)||^,_3(„) <i^At, te [o,r], 

where K depends on s, T and uq only. Here v(t, t) denotes the splitting approxi- 
mation in our approach evaluated on the "diagonal" (t,T) = (t, t). 

To obtain second order convergence it is common to apply the Strang splitting 
formula, thus formally 

($B(At/2)$^(At/2))($^(At/2)$B(At/2))j =0{Af), 

where t„ — > t as At ^ 0. Here we show rigorously this result for the KdV equation. 
Indeed, in Theorem |3.5| we prove that if uq G for some s > 17, then for At 

sufficiently small 

\\v{t,t)-u{t)\\H^-o <KAt'', tG [o,r]. 
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where the constant K depends on uq, s and T only. Again v(t,t) denotes the 
operator splitting approximation in our approach. Observe that we have to increase 
the regularity of the initial data, and hence of the solution, in order to get increased 
accuracy of operator splitting. Note also that with this type of operator splitting 
our error estimates are in a much weaker norm than the assumptions. 

It is clear that the present approach applies to several other equations, and this 
is currently being investigated. Furthermore, for applications to numerical analysis, 
one would need to replace the exact solution operators and by numerical 
approximations, say <i>^ and <i>^, and study their behavior in the limit as (5 — > 0, 
and also to replace the time derivatives by discrete differences. Again this will be 
studied separately. 

The paper is organized as follows. We start by presenting the Godunov operator 
splitting method for abstract operators. Next we apply this approach to nonlinear 
ordinary differential equations where the procedure is fairly transparent, before 
we discuss Godunov splitting for the KdV equation. Subsequently we present the 
Strang splitting for abstract operators. Finally we apply this procedure to the KdV 
equation. 



We first present a formal calculation motivating the rigorous analysis that will 
follow. Consider an abstract differential equation 



for some fixed positive time T, where C typically will be a differential operator in 
the spatial variable. We assume that uq and u are in some Hilbert space X, and 
write the solution as 

u{t) = <i>cit;uo). 
Formally, expanding the solution in a Taylor series we find 

u{t) = uo + tutiO) + ©(i^) = uo + tC{uo) + 0{f), 
where we have used the equation Assume that one can write 

C = A + B 

in some natural way. Operator splitting (of the Godunov type) then works as 
follows: Instead of solving the problem with C directly, one alternately solves for 
small time steps the equations Ut = A(u) and Ut = B(u). Making the time steps 
finer and finer, the approximation will presumably converge to the solution of the 
original equation Q. 

More precisely, fix a positive time step Ai, let i„ — nAt, n e No, and define a 
family {uAt(tn)}n of functions 

UAt{tn+i) - '^A{At; $B(At; UAt(i«))) = ^A{At) o $B(At)MAt(t„), n e No, 

UAt(O) = uo- 

The traditional method of extending the solution to any t e [0,r] has been to let 
"time run twice as fast" in each of the subintervals [tn,tn+i/2\ *md [tn+i/2i ^n+i], 
where i„+i/2 — tn + At/2. Thus obtaining 



2. Operator splitting 



(2) 



U|t=0 = Uq 
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in some norm. 

The convergence can be improved to second order by using the Strang sphtting 
formula. To this end we let the approximation (this time denoted v^t to distinguish 
from the previous approximation) read 

) - $B(At/2; <^A{At; ^B{At/2; VAtitn)))) 

= $B(At/2) o $^(Ai) o $B(At/2)wAt(i„) 

= ($B(AV2)o*A(Ai/2)) 

o ($^(Ai/2) o $B(AV2))t^At(i„), n e Nq, 

UAt(O) = Uo- 

Formally one now has 

\\vAtitn) - u{tn)\\ < O(Ai^) as Ai ^ and t„ -> t. 

To show that the operator splitting solutions are well defined, we shall later make 
use of the following bootstrap lemma, taken from [8 , Prop. 1.21]. 

Lemma 2.1. Let t G [0,1"]. Consider a continuous function (j): [0,T] — > [0,cxd). // 
there exists a positive constant a such that: 

(a) m < a; 

(b) for any t such that (j)(t) < a, we can show that (f)(t) < a/2; 
then 4){t) < a/2 for all t e [0,r]. 

2.1. Doubling the time variable. We shall formulate the operator splitting so- 
lution by introducing two time variables, and define a function v = vit, r) for it, r) 
in the set 

LT/AtJ 

^At = IJ [tn^tn+l] X [t 

by requiring that 

i;(0,0) = uo, 

(4) Vt{t,tr,) = B{v{t,tr,)), te{t„,t„+i], 

Vr{t,T) = A{v{t,T)), (t,T) G [t„,tn+l] X (t„,t„+i], 

where n = 0, . . . , [T/AiJ . Observe that 

UAt{tn) = w(t„,t„), n = 0, . . . , \ T/M\. 

where ma* is given by (|3|. The specific extension v of {wAt(^n)}n to [0, T] will serve 
as an important technical tool to be utilized in the analysis. 

The exact solution of ^ is still denoted by u. Introduce the error function 

w{t) = v{t,t) - u{t). 

The aim is to show that 

\\w{t)\\<0{M), At^O, te[0,T] 

in an appropriate norm. 

We introduce second order Taylor expansions of the operators A and B (see [TJ 
p. 29]), 

A{f + g)=.A{f) + dA{f)[g]+ [ {1 - a)d('^ A{f + ag)[gfda, 

Jo 

B{f + g) = B{f)+dB{f)[g]+ [ {1 - a)d^'^ B{f + ag)[g]^da. 

Jo 
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Figure 1. Schematic view of the Godunov sphtting and the defi- 
nition of v{t,T), cf. equation Q. 



Thus we find 

Wt — dA(u)[w] — dB(u)[w] — vt + v^- — Ut — dA(u)[w] — dB{u)[w] 

= vt + A{v) - {A + B){u) - dA{u)[w] - dB{u)[w] 
= vt- B{v) + {A{v) - A{u) ~ dA{u)[w]) 
(5) + {B{v) ~ B{u) ~ dB{u)[w]) 

= F{t)+ {l-a)d'^^'^A{u + aw)[w]'^da 



(1 - B{u + a'w)[w]'^ da 

where we have introduced a forcing term F{t) = F{t, t) defined as 

(6) Fit,T)^Vt{t,T)-Biv{t,T)). 

We can rewrite (|5| as follows 

(7) Wt-dC{u)[w]^F{t)+ [ {l-a)d^^^C{u + aw)[w]^da. 

Jo 

The forcing term satisfies the following time development 
F, - dA{v)[F] = Vtr - B{v)r - dA{v)[vt - B{v)] 

= A{v)t - dB{v)[vr\ - dA{v)[vt] + dA{v)[B{v)] 
= dA{v)[vt] - dB{v)[A{v)] - dA{v)[vt] + dA{v)[B{v)] 
= [A,B]{v,v), 
where we have defined the commutator 

[A,B]if,g) = dAif)[Big)]^dB{f)[A{g)]. 
For simplicity we will subsequently be writing _B](u) rather than [A,B]{v,v). 

2.2. Ordinary differential equations. As a warm-up we consider the ordinary 
differential equation 



(8) 



(9) 



ut = C{u), t > 0, m(0) = uq e 
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To simplify the presentation we assume that C is a quadratic function, i.e., that 
C" = cPC = or (PC[f, g] is constant. This means that the integral in ([t]) reduces 
to a constant. Furthermore, we assume that (|9| is such that there is a unique 
solution u{t) such that \u{t)\ < K^q^t for t G [0,T]. Furthermore, we will assume 
that the operators A and B are two times continuously differentiable and 

(10) A\A"eL°°, A(0)=0, \B{u)\<K\u\'' , \B' {u)\ < K \u\ . 

Throughout this paper we use the convention that for a quantity a, Ka denotes 
a constant depending on a (and perhaps other things). We use this notation to 
highlight the dependence on a. The actual value of Ka may be different at each 
occurrence. 

Let now a be a positive constant (its precise value will be fixed later) . To start 
the bootstrap argument we assume 

(11) \v{t,T)\<a, {t,T)enAf 

In this example, since C is quadratic, (PC — k for some constant symmetric matrix 
K, and w satisfies 

(12) wt + C'{u)w = F + w'^'^w, t>0, w{0) = 0. 
Furthermore, F satisfies 

Fr - A'iv)F = [A, B]{v), {t, t) e [i„, i„+i] X (i„, i„+i], 
and F{t,tn) = 0, for each n. This means that 



d 



\F\ <K\F\+Ka 



2 



dr 

Hence, Gronwall's inequality implies that 

\F\{t,T) <KaAt, 



where Ka is a constant depending on the assumed bound on u in ( 1 1 1 . In view of 

w\ < Ka \w\ + KaAt, w{0) = 0. 



this bound and (12 I, 

d 
dt 



Gronwall's inequality gives that 

\w{t)\ < e^'-HKaAt < KaAt. 

Trivially we have 

\v{t,t)^v{t,T)\< \A{v{t,s))\ds<KaAt, 

J min{f,T} 

for (i, r) e [ir!,in+i] ^ {tmtn+i] for Siiiy n. Then we can conclude that 

\v{t,T)\ < \u{t)\ + \w{t)\ + \v{t,t) - v{t,T)\ <K+KaAt. 

Now we are in a position to choose a so that K < a/4, this determines Ka- Next 
choose At so small that KaAt < a/4, then 

\v{t,T)\<a/2 (t,r)Gf]At. 

Hence, by the bootstrap lemma, |u(t,t)| < a/2 for all t. Consequently, 

(13) \v{t,t) - u{t)\ < Ka/2At, 

i.e., the operator splitting is as expected first order accurate. 
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EiTOr (u0=.9, N=30, At=.l) 




Figure 2. (Left) The exact (blue) and the approxuiiate (red) so- 
lution. (Right) The error. (Below) The error w{l) as a function of 
At. 



Remark 2.2. An interesting example is the logistic equation u' = u(u — 1) where 
we can write A{u) = —u and B(u) = . Exact solutions are available for all 
operators involved, specifically 

*c(t)"o = —l^. r, ^A{t)uo = uoe-\ <i>B{t)uo = ^ ^° . 

Uq + e*(l - uq) 1 - uot 

Let uq (z (0,1). Then there is no blow-up in the full equation, but blow-up for the 
equation Ut = b(u) at t* = I/uq- The function v reads in our case 

^(^'")=l-.(t„,t„)(t-t„)' 

where 



(1 - e-^*)e*" + MoAi(l - e*-) ' 

The quantity v{tn,tn) is well-defined on [0,T] if one chooses At such that 

1 - e~^' 

At < 



uo(l-e-*")' 

Since t^ < T, and < Uq/(uq — 1), we find that uq{1 — e~*") < uo(l — e"-^) < 1. 
Thus we have to choose At such that 

.o(l-e--)<^(l-e-^*) = l-^ + 0(At^), 

or 

Ai< 2(1 -1*0(1-6-^)). 
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In this case one can verify (13) directly, namely 

\v{tn,tn) - u(<„)| 

^1 |2.t„_,^ \m'\l~e-^^)-l\ 

^|K + e*"((At)-i(l-e-^*)-"o))K + e*"(l-«o))| 

= ©(At). 

T/ie example is illustrated in Figure^ 

2.3. The KdV equation. Let us now apply this general framework to the KdV 
equation, that is, 

ut = uu^ - Uxxx, u\t=o = Mo e H'^iR). 

In this case C(u) = uwj; — Ur^^xi f^nd the evolution operator ^c{tj ' ) ■ H"^ H'^ is 
bounded. Therefore, the Hilbert space M" of the previous example is replaced by 
the Sobolev space _ff*(IR) with the inner product and norm 

(/,5)h-'(R) = / dif{x)dig{x)dx, = (/, 

We choose A to equal (minus) the Airy operator: 

^(y) fxxx: 

dA{f)[9] = -9xxx, 
d^'U{f)[g,h] = 0, 
and B to equal the Burgers operator: 

Bif) - //., 
dB{f)[g] = fg, + f,g, 

d^^^B{f)[g,h]^hgx + hxg, 
S'^B{f)[g,h,k] = 0. 
In this case the commutator reads 

[A,B]{fJ)^-^dl{fxr. 
Thus the equations ([7|, ([6]), and ([s]) are 

(14) Wt - iuw)x + Wxxx = F + WWx, 

F = vt~ vv^, 

(15) F, + Fxxx ^ ~ldUvx)^ 
respectively. 

From we recall the classical result: For uq G H''{R) with s > 2 there exists a 
unique solution ?i e C([0, T], //''(M)) of the KdV equation 

In particular, we can assume that there exists a constant K (depending on T, uq, 
and s) such that 

l|w(t)ll^f.(K) <i^, te[o,r]. 

To save space and typing efforts, we shall write H'^ for H^{M.), and d for 9a;. 

Next, let s be an odd integer larger or equal to 5. Assume that there exists a 
constant a such that 

\\v{t,T)\\^i, <a, (i,T)ef]At, 

where k{s) = (s — l)/2. 
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Let us estimate the behavior of the Airy and Burgers operators. The Airy 
equation leaves all Sobolev norms invariant, viz. 

Mt,T)\\H^^ = l|w(i,in)lliffe ■ 

By definition we find for the Burgers operator 
2di 



n — n I — n \ / 



j=0 k=0 

For j < s, any of the above terms can be estimated by 



dx 



^max{fe+l j-fe}^ ^min{fe+l j-fc}^ 

L2 



L2 



<K\\vrH^,\\v\\j^, 



<Kc.\\v\\h.. 

For j = s all terms with k < s can be estimated similarly, 



< Wd'vW 



L2 



^inin{fc+l,s— fc}^ ^max{fc+l,s — fc}^ 



<K\\v\\%^,\\v\\^, 

<K^\\v\\l.. 
We are left with the term where k = s — j, viz 

d'vd'+^vvdx ^ 



{d%y dvdx 



<K\\v\\^, Ml. 



L2 



if A; > 2, i.e., if s > 5. 
Thus 

(16) 



dt 



\\v{t,tn)\\Hs <K^\\v{t,tn)\\Hs, 



which implies that 

In particular, for any n, 

\\vitn,tn)\\H^, <e^°^*|lt;(i„_i,<„_i)|l^, < e^°*" Ikollff. 
Thus we have shown the following result. 
Lemma 2.3. We have 

\\v{t,T)\\^. <K^, {t,T)enAf 

Observe the general result, obtained by integration by parts. 



Next we analyze the forcing term that satisfying (15 1. By taking the inner 
product with F in ( 15 1 we get 

\dr \\F{t, r)|l^._3 = -(F, F^^^H'-. - ^(F, d\v,f)H^^^ 
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Since H'^ is an algebra 

ll^'^'(^-)'L=-3 < \\{v.f\\Hs-^<K\\v£H^-^<Kl 

Also F{t,tn) = for t e [tn,tn+i]; and we conclude that 

\\F{t,T)\\^,_s<K^At, (t,r)er!At. 

As for the error function w — v{t, t) — u(t), we have the following estimates. Let 
E{t) = \\w{t)\\^s-3- By taking the inner product with w in ([l4|) 



1 d 
2dt 



(17) 



= (w, {uw)x)h—^ + {w, wWx)h---^ + (w, F)jis-3 

s — 3 „ 



The first integrand on the right is expanded by Leibniz' rule. We get 

mi ( ) / d'wd''+^wd^-''u + &wd''w&+^-^u + d^wd''wd'+^-''w dx. 

j=0 k=0 ^ 

For 0<fc<j<s — 3, we can estimate 



<K\\wr^^,^s\\u\ 



w\ 



For J = s — 3, we can use the same strategy for those terms with fewer than s — 2 
derivatives on w. The term with s — 2 derivatives on w can be estimated as 

^ {d'-^wf dudx < Xllull^, ||u>||^,_3 , 



9" ^wd'^ ^wudx 



(9" ^w) dwdx 



The last term in ([T7| is overestimated by ||F||^s_3 ||u)||^s_3, and we get 



dt 

which implies 



E'^{t) < K^E\t) + KE{t) ||F||^,_3 < K^E'^it) + KM, 



dt 



E{t) < KaE{t) + Kilt. 



Since -^(O) — 0, Gronwall's inequality yields 

\\w{t)\\H^^a =E{t) < Ka^At. 
Recall that t)||^^_3 = t)||^s_3 because 

9t 11^.-3 = -2{v,Vxxx)hs = 0- 

Now u is bounded in i?", and we infer that 

||v(i,T)||^,_3 <K + K^At. 
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Since s > 5, 

s - 3 > fc = 

and we get 

\\v{t,T)\\^t <K + K^At. 

First choose a > AK, then choose Ai so small that K^At < a/4, so that K + 
KaAt < a/A + a/4 = a/2. Hence, by the bootstrap argument we have proved the 
following theorem: 

Theorem 2.4. Fix T > 0. Let uq e H'^iR) with s > 5. Then for At sufficiently 
small we have 

||«(i,t)-w(t)||^._3(R) <ifAi, [0,T], 

where K depends on s, T and uq only. 

Remark 2.5. Instead of defining v by ^ for the KdV equation, we could also 
interchange the order of the Airy operator A and the Burgers operator B in the 
definition of v. The same procedure as described above would apply, and Theorem 
\2.4\ would remain valid. This remark is important for the Strang splitting to be 
discussed next. 

3. Strang splitting 

To achieve higher-order convergence it is common to consider the so-called Strang 
splitting. Now we approximate the solution by using two Godunov splittings, each 
with a time step of At/2, and in alternating order. Explicitly, we define 

i;(0,0) = Uq, 
Vf{t,tn) = B{v{t,tn)), t e (t„,t„+l/2], 

(18) Vr{t,T) = A{v{t,T)), {t,T) € [i„,t„+i/2] X (t„,t„+i/2], 
1'r(t«+l/2,T) = ^(w(in+l/2,T)), T G (i„+l/2, ^n+l] , 

Vt{t,T) = B{v{t,T)), {t,T) e {tn+l/2,tn+l] X [i„+l/2, ^ri+l] , 

for n — 0, . . . , [T/AtJ . We consider this function for [t, r) in the domain 

[T/At] ^ 

f^At = [J i[tn,tn+l/2] U [t„+i/2 , tn+l] j- 
ri=0 

The aim is now to show that 

\\v{t,t) -u{t)\\ =0{At^) 

in an appropriate norm. Here the abstract analysis of Section |2.1| applies, and we 
find: 

(i): On the domains [in,in+i/2] x [tn,tn+i/2] we have 

Wt-dC{u)[w]^ F{t)+ [ {l-a)d^^^C{u + aw)[w]^da, 
Jo 

(19) F{t,T)=vt{t,r)-Biv{t,T)). 

The forcing term satisfies the following time development 
Fr^dA{v)[F] + [A,B]{v), 
Ft = Vtt - dB{v){F\ - dB{v){B(y)\. 
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Figure 3. A schematic view of Strang splitting. 



(20) 



(ii): On the domains [tn+i/2,tn+i] x [tn+i/2jtn+i] we have (here we write Wt 
rather than Wr since m; is a function of one variable only) 

Wt-dC{u)[w] = G{t)+ f {l-a)d^'^^C{u + aw)[w]^da, 
Jo 

G{t,T)=Vr{t,T)-A{v{t,T)). 

The forcing term satisfies the following time development 

Gt^dBiv)[G] + [B,A]{v), 
Gr = vrr - dA{v)[G] - dA{v)[A{v)]. 



We extend F and G to all of f^At using the same definitions, (19 1 and (20 1, respec- 
tively. Observe that this implies that = on [tn+i/2,tn+i] x [tn+i/2,tn+i], while 
G = on [tn,tn+i/2] X [^ri , ^n+i/2] ■ The total forcing term is defined by 

H{t,T) = Fit,T)+Git,T). 



3.1. Ordinary differential equations. One can consider the case of ordinary 
differential equations, as wc did in Subsection |2.2| for Godunov splitting, but for 
reasons of brevity we will only revisit the example in Remark |2.2[ 



Remark 3.1. We find, using the definitions (18), that 

''^(Wl/2,«„ + I/.) fort TG\t , f 1 

for n ~ 0, . . . , \ T/At\ . By induction we determine 



(22) w(t„,i„) 



(1 - e-'^*)e*" +uoAi(e*" - l)(e'^* + l)/2' 



OPERATOR SPLITTING FOR THE KDV EQUATION 
u'=u(u-l), u0=.9 



EiTor(uO=.9, N=30, At=.l) 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 



0.5 1.0 1.5 2.0 



Figure 4. (Left) The exact (blue) and the approximate (red) so- 
lution. (Right) The error v{t,t) — u{t). 



In this case we compute, when we for convenience write a = Uq(1 — e*"), 

\w{tn)\ = \v{tn,tn) ~ u{tn)\ 



uo(l-e-^') 



Uq 



(23) 



(1 - e-'^*)e*" + woAi(e*" - l)(e'^* + l)/2 e*- + uo{l - e*-) 
1 1 

\uo\ 



< 



e*" + aAt{l + e^*)/(2(l - e^^*)) e*- + a 
|uo|'e*"(e*"-l) 



|(e*" + aAt{l + e^*)/(2(l - e-'^*)))(e*" + a)\ 



At{l + e 



At^ 



2(1 - e-^*) 



The example is illustrated in Figure^ Observe the strong oscillations in the error, 
it is these oscillations which prevent the error from growing too large. 

3.2. The KdV equation. For the KdV equation, we use B{v) = vvx and A{y) = 
—Vxxx, and the above analysis yields: 

(i): On domains [i„,i„+i/2] x [i„,i„+i/2] we have 

Wt - {uw)x + Wxxx = F{t) + WWx, 
F{t,T) = Wt - VVx- 

The forcing term satisfies the following time development 

3, 

= -tl^^ - 

(25) 

Ft = vtt - {vF)x - {2vvl + v^Vxx), 
since dB{v)[B(v)] — v(vvx)x + Vx{vvx) — ^vv"^ + v^Vx 



(24) 



(26) 



(27) 



(ii): On domains [t„+i/2, in+i] x [t„+i/2, we have 

Wt - {uw)x + Wxxx = G{t) + WWx, 

G{t,T) = Vr{t,T) +Vxxx- 

The forcing term satisfies the following time development 
Gt = ivG)x+ldUvxr, 

G-j- ^rr ~l~ ^xxx V^- \ 

since dA{v)[A{v)] = (vxxx) XXX ■ 
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To start the bootstrapping procedure we fix an odd integer s and a positive 
constant a, whose values will be determined in the course of the argument. Now 
assume that 

where fc=(s — l)/2. As a consequence, we have 



cf . the proof of Lemma |2.3| which can be easily adapted to Strang splitting. 
We need to introduce the function 



(28) 



z{t) = w(t) + w{t +^)= w(t) + w{t). 



In the following we will write = (j){t + ^) for any function (j). Straightforward 
calculations yield that z satisfies 

,1 



(29) 



zt ~ {^z^ + UZ - Zxx)x = H + H + {w{u -u) - ww)^^ 



By our techniques, we must work in ^ (so at least s > 9), therefore set 

E{t) = Mt)\\Hs-.i^n)- 
By taking the i?*^^ inner product with z in (p9|) we get 



^^-B^(<) = [z, i^z'^ + UZ- Zxx)x)h^-o + {z,H + H)hs-9 + (z, {w{u ~ u) - ww)x)h 
— (2, zzx + UxZ + uZx)h=-^ + {z,H + H)ff!,-o — (zx,w{u ~ u) — Ww)hs^9 



< (Z, ZZx + UxZ + UZx)h=-B + E{t) 



+ E{t) \\w{u — u) — ww\\ 



< K^E^{t) + KE{t){ 



H + H 



H + H 

+ \\w{u — u) — WwW^s-g ) 



Thus by Gronwall's inequality 

ft 

E{t) < E{0) 



( H{a,c7) + H{a,a) 



(30) 



+ \\w{a){u{a) - u{a)) - w{a)w{cr)\\jj,-o ) da 



< E{0) 



H{a, a)+H{a,(7) 



+ \\w{a){u{(T) — u{(t)) — w{a)w{a)\\f^s-g ) da. 



Next we turn to the detailed estimate of each of the terms in ( |30| . We start 
with most involved one, the forcing term, which can be estimated as follows. We 
consider the term F first. Since F{t,tn) = 0, we easily see that Ft{tn,tn) = 0; thus 

(cf. ^) {Ft + Fr){tn,tn) 



id^M^{tn,tn). Thus 



F{t,t) = {Vt ~ VV^){t,t) 

= F{tn, U) + {Ft + Fr){t„,t„){t - t„) 

, {t-tn)^ 



{Ftt + 2Ftr + Frr){cr{t ~ tn) + i„, a{t ~ t„) + tn)da 



= -l^dl{v,f{tn,tn){t-tn) 



{t-t^f 



{Ftt + 2i^tr + Frr){(T{t - U,) + t„,CT(t - t„) + tn)da. 
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As for the second derivatives, we find 

3 

^dlF+\dl{v,f-?,dl{v,v^''^), 
{61) ^ 

■^rt '^ttxxx {'^xxx'^x ~l~ '^'^xxxx^ti 

Ftt = vttt - VttVx - 2vtVxt - vvxtt- 

Similarly we find for the forcing term G the following estimates. Since G'(t„+i/2, t) = 
0, we easily see that G'^(t„+i/2,i„+i/2) = 0; thus (cf. ([27])) (G't+G^)(i„+i/2, t„+i/2) = 

ldl{Vxf{tn+l/2,tn+l/2)- ThuS 

At At, , A< At, 

G{t+—,t+—) = {Vr+ V.,,,){t + — , i + — ) 

— G(t„+i/2; in+1/2) + {Gt + Gr){tn+l/2, in+l/2)(^ — tn) 
(t-t )2 /-l 

+ ^ ^ / {Gtt + 2Gtr + Grr){ait-tn)+tn,ait-t„)+tn)da 

^ Jo 

3 

(t-t V 

+ ^ / {Gtt + 2Gtr + Grr)Ht~t,,)+t„,a{t-t^)+t^)da. 

The second derivatives Gu + ^Gtr + Grr will have to be considered similarly to 
those for F. These read 

Gu = 2vlG + ivu^G^ + v^G,:^ + Gvv^^ 

3 

Grt = VrrVx + 2VrVrx + VVrrx + 9'' (Gi;) + -8^ {yl) - 8^ (Wx) , 

GfT ^TTT ^" ^TTXXX- 

Lemma 3.2. We have the estimate 

mt,T)\\l,_, + \\^{t,f)\\l^_, < K^, (i,T),(i,T) el^At, 

where $ = {Fu + 2FtT + Frr) and * = {Gu + 2Gtr + Grr)- 



Proof. We have that ||F||j:^s_9 < K^At, we shall get a similar estimate for |lG||^s 
For tn < T < i„+i/2, we have that 

Gt^VxG + vG^ + ^d^ (vl), 
and G{t,tn) — 0. Taking the iJ*^^ inner product with G we get 

s-9 



7 n J 



fc=0 ' 

The first term expands by the Leibniz rule; a typical term in this expansion reads 
(here < j < A: < s - 9), 

j d^+^vd^-^Gd'^Gdx < ||<9^'+H'||i^ l|G||^s-9 . 

Similarly the second term can be expanded and estimated, except the term 
containing d^^^G, which is estimated as 

J vd^+^Gd^Gdx {d^Gf dx = j dv {d^Gf dx 
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<l\\dvh^ ||G|Ih.-9. 



We bound the last term as 

d'^+^vl)d''Gdx <i^||^;||^.||G||^,_.. 



Summing up, we get 

1 1|G||^3-. < \\G\\l.-. + \\G\\^s-. . 
Using Gronwall's inequality and that G{tn,T) = we get 

< K^At. 

We also need estimates for vt, Vu and vai. where t,r G [in,in+i/2]- In this set, 
Vt = —Vxxxi and this evolution preserves the H'^ norm. For r = 

Vt = VVx, 
Vtt = VtVx + VVtx, 
Vm = VttVx + 'iVtVtx + VVttx- 

Therefore 

\\vt\\H<^ < WvWfjk \\v\\Hk+i < K\\v\\%k+i , 

WviuWh'' ^ ^ II^IIh'=+3 ■ 

Next we turn to estimates of Vr, Vtt and Vttt in the set i,T G [ira+i/2) in+i]- Here 
Vt = vvx, and setting ui = Vr, = Vrr and 7 = Vrrr, we get 

Wt = VxU) + vux, w(i„+i/2, r) = -d^v, 

6t = 2uJUIx + VxO + v6x, 0{tn+l/2, t) = d^v, 

7t = 3eux + BOxOJ + Vx'y + vjx, 7(*n+i/2, r) = -d^v. 

Starting with w, for A; < s — 9 we get 

fc 



ld_ 
2di 



dx. 



Using Leibniz' rule, all terms except the ultimate one will be of the type 



/ 



< \\v\ 



The ultimate term (with one too many derivatives on u) is as usual estimated as 



d^+^ujd^ujdx 



dx 



< \\v\ 



II l|2 



Gronwall's inequality then yields 

Reasoning similarly for 6, we find that 

II^IIh. it) < e^-^* ( max ||u;||^.+, (s) + \\0\\j,, (1^+,^,] 

< {\\v\\hI'+^ + Mnk+e^ {t„+l/2,r) . 

Finally, the estimate for 7 reads 

+ 11^11 iJfc+g) {tn+l/2,T) . 



OPERATOR SPLITTING FOR THE KDV EQUATION 



Summing up 
Now 

Working in the square [in, t„+i/2]^, from (31 1, 



d^dvf +3\\d'^ {dvd^v)\\j^^ 



(32) 
and 



\\Ftr\\H^-9 < \\vtt\\H'^-<^ + Ka{\\v\\jj^,-s HvtH^.-e + ||wt||^.-9 \\v\\hs-6 
(33) <^a(||«||^.-. 



V\\hs-7 \\V\ 



and also 

(34) < K^. 

Hence ||$||j:^,_9 < K^. 

In the second square [iTi+i/2j ^n]^, we write 

II^IIh^-^ < \\Gtt\\H^-9 + 2 |lGt.||^._o + \\Grr\\Hs-9 ■ 

Each term above is estimated individually as 

K^{\\v\\h.^s ||G||^._9 + Mh^-s ||G||^._« + IIHIh.-o ||G||^. 



(35) 



||Gt^||^,_9 < ||Wrr«x|lif.-9 + 2 |lWrWx2;|!H=-9 + II ««rrx || .-9 + 1 1 9'* (Gw) 1 1 

+ ^ll^'(^')L=- + l|5'(««-)L=-« 



<Ko,(\\v 



+ I|G||h.-5 ||«||h.-5 + ||«|Ih.-3 + ||v||h.-3 ||«||h.-2 

(36) < K^, 

IIG-rTll/fs-g 5: || '^'ttt || s-9 + || "yrTa^ica; || //s-9 
<^c,(||«|lH=+ll"llff=) 



V H-.-9 V 
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(37) < K^. 



Collecting ( 34 1 , ( 33 1 , ( 32 1 , (351, ( 36 1 and ( 37 1 finishes the proof of the lemma. □ 



Lemma 3.3. The map 

is Lipschitz continuous with Lipschitz constant at most K^. 
Proof. Set w{t,T) = d'^{vl){t,T). Then we have 

'^t — '^^xx^txx ~t~ '^Vtxxx'^x ~t~ '^Vxxx'^tx ^ 
'^T — '^'^XX^TXX ^" '^^TXXX^X ^" '^^xxx'^Tx • 

In the square [tmtn+1/2]'^ we have Vr — —Vxxx, thus 

Wr = —'iVxxd^V — 2d^VVx — 2Vxxxd*V. 

Hence in this square 

\\wt+Wr\\Hs-9 < K(^\\v\\jj,^r \\vt\\H^--' + Iktllff-s \\v\\hs-s 

+ II^^IIh=-6 \\vt\\H^-s + MH^-r \Mh^-^ 
+ + ||l'||^.-6 ||w||^,-5 

In the second square [tn+i/2^tn+i]'^ we have Wt = vvx, and 

Wt = ^Vxx {vVx)xx + 2 {'VVx)xxx '"^ + '^Vxxx {vVx)x ■ 

Therefore in this square 

\\wt + Wr\\Hs-9 < K(^\\v\\fjs-7 \\vVx\\h^-7 + \\vVx\\hs-6 IK'Uffa-g 

+ \Mh^-<^ W^VxWhs-s + \\v\\hs-s \\Vt\\h^-» 
+ Ikllff-s ll^^r|lH=-6 + \\v\\hs-6 ||Wr||ff.-8^ 

< ^(ll^IlL-- Ikllif-e + MIh'^--' \Mh'^-<^ \Mh^-s 
+ \Mh^-o \Mh^-^ IMh^--' + Ikllff-s 
+ \Mh^-<^ ll«llif-3 + \\v\\H,-e \\v\\h,- 



□ 



From this lemma it follows that 

m = -d\vxf{t,t)+d\vxf {t+Y^'+Y 

satisfies 

|ir(i)|l^._, <i^„Ai, te[o,T]. 

The following lemma will be convenient: 
Lemma 3.4. For t G [tm,tm+i/2\ we have 

(38) / (^^ (^W + G{t + At/2)))^ dx < K^M^ 

3=0 

(39) / (^^ (^W + G{t- At/2)))^ dx < K^At^ 
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Proof. Wc show (38 1, (39) is proved similarly. By a Taylor expansion, for t € 

[tm, tm+1/2] 
s-9 „ 

/ {d^F{t) + G{t + At/2)))^ dx 



s-9 



2 70 



$ (CT(t - tm) + t,n) da 



tm ) 



2 (-1 



* {a{t - tm) + tm) da^ 



dx 



3=0'' 



{t tm) 



2 ^1 



(a(i - tm) +tm) + -^ {cr{t - tm) + tm)) rfcr) 



dx 



-Ay.! (-^' (^-)' + (^-)' + ^^2)) it ~ tm)y dx 

3=0-' 

+ 2^y" (^d' {<i>{a{t~tm)+tm)+'^{<j{t-tm)+tm))dait-tm)^^ dx 



(^'^ i<i>Ht-tm)+tm)+'^i<jit^tm)+tm))da^ dx {t - tm)' 



< K^At^ {t-tmf 



„1 

+ 4 / (\\^{ait-tm)+tm)\\ls-o + \\'^ia{t-tm)+tm)fHs-9)da{t-tm)' 



< K^At^ [t - tm) +K^{t- tm) 

< K„At^. 



Combining the above results, we find that for t £ 



□ 



^ ||ff(CT) + H((T)||^^_^ d(T < ^ " H{a) + H{u) 



Hs-9 



da 



j-t„, + i / "-"^ r/ A. x2 ^l/2 



m=0 * 



m=0 * 

n— 1 

<^E 



s-9 



J=0 



At 



At. 



E/ E / + + dx] da 



1/2 



s-9 



+ (a.(F(^+^) + G(a))) 



2 \l/2 

da; I da 
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< 



m=0 ' 



1/2- 



.1 + 1/2 



s-9 



1/2 



1 + 3/2 



^ / {d^ (F(cr) + G((7 + At/2))) dx \ da 

s-9 „ ^ 1/2 

/ (a^' (F((t) + G((T - AV2)))' dx ) da 



< i^aAi^ E / + E 



da 



where we have used Lemma |3.4[ This finishes the estimate for the forcing term. 

Next we turn to the estimate of the term J^* || w(cr)(it(tT) — — w(a)w{a)\\ j^s-'^ da 

in ( 30 ) . Here we can use the estimates from the Godunov sphtting to infer that 

ll«^(CT)llff.-9 + < ^aAi, <y e [0,T]. 

From the KdV equation we infer immediately 

/•At/2 

||u((t) - u((t)||^,_9 < / \\(uux- Uxxx){o + T)\\^,-^dT <Kbd, cre[0,T]. 

^0 

Thus 



\w(a)(u(a) — u{a)) — w((7)ti;(cr)||^s_9 da 



< K 



(||w(cr)||^,_9 ||(u-u)(cr)||^,_9 + |1w(ct)|1^,-9 \\w{(^)\\h^ 



da 



(40) < Kc.M'^ 



The last term to estimate in ([30]) is £:(0) = w(0) + w{Q) = w(A</2). For 
i < At/2, we find that 

w(t) = i) - u{t) = t [ {B{v{st, 0)) + A{v{t, st)) - C{u{st)))ds 
Jo 

<.i 

t 



(^B{u{st)) + J dB{u{st)+a{v{st,0)-u{st)))[v{st,0)-u{st)]da 

+ A{u{st))+ I dA{u{st) + a{v{t,st) - u{st)))[v{t,st) - u{st)]da 
Jo 

- {A + B){u{st))'jds 

(^dB{u{st) +a{v{st,0) - u{st)))[v{st,0) - u{st)] 
+ dA{u{st) + a{v{t, st) - u{st)))[v{t, st) - u{st)]) ds da 



1 pi 



= t 



^0 



dB{u{st) + a{v{st, 0) - u{st))) 



dr 



iv{T,0) - u{T))dT 



f d f d 1 \ 

+ dA{u{st) + a{v{t,st)-u{st))) / —{v{t,T)-u{T))dT+ — i;(t,0)o?t) 

I In dr ./n dr \l 



ds da 
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dB{u{st) + a{v{st, 0) - u{st)))[B{v{T, 0)) - (A + B){u{t))] 



+ dA{u{st) + a{v{t, St) - u{st)))[B{v{t, t)) ~{A + B){u{t))] ) dr 



+ I dA{u{st) + (7{v{t,st) ~u{st)))[B{v{T,Q))]dT 



= t 





r-l i-l 

In Jo 



ds da 



{u{st) + a{v{st, 0) - u{st))) {v{t, 0)v{t, 0)^ - u{t)u{t)^ + u{t)^^^)) 



+ / Ut,0)v{t,0)^ 

In \ / XXX 



)dT 



ds da. 



Taking the H^^^ norm above, using the triangle inequahty, the bounds on v and u 
and the Cauchy-Schwarz inequahty, we find that each of the above integrands are 
bounded by Ka- Thus 



(41) 



\\wmH^-o<t 



/ K^dT+ / 

'-JQ Jo 



lO Jo 



Kndr 



ds da 



Hence we infer that 

(42) E{Q) = \\w (At/2)||^,_3 < X„At2 



Collecting the estimates from (40 1, (41), and (42 1, we find that (30) reads 

(43) E{t)^\\z{t)\\u,^,<K^At'. 
By the triangle inequality 

(44) \\w{t)\\„.-, < 2 |lu;(t)|l^.-o < lk(i)llff.-o + \\w{t) - w{t)\\^,^, . 
To estimate the last term on the right-hand side we write: 

\\w{t) - w(i)||^,_9 < ||w(t„) - w(i„+l/2)||^.-9 + \\w{t) - w{tn)\\Hs-9 

+ ||u;(i + At/2)-u;(t„+i/2)||^_9, 
for t e [t„, tn+i/2]- (Similar expressions hold when t G [tn+1/2, tn+i]-) We note that 



(45) 



(46) 



* ,At, .At. , , ^ .At. . . 

w(t„+l/2) = $A(^)«'B(^)w(i«,i«) - $c(^)M(in), 

W{t) = <^>A{t - tn)<^B{t - tn)v{tn,tn) - ^cit - tn)u{tn) , 

W{t+^) = «'A(i "irO*B(<-^n)^'(in+l/2,in+l/2) "^c(^)u(i„+l/2), 



when t e [tn,tn+i/2\- Each of the expressions on the right-hand side of (45) needs 
to be estimated: 

wit.n+1/2) - w{t„) = <^A{^;'^Bi^;v{tn,tn))) - <^ c ] u{t„)) - (w(t„,i„) - u(i„)) 

^ .At ^ .At . ... ^ ,At . 
= $A(^;«'s(^;w(i«,irO)) -*c(^;i'(i«,in)) 

+ ^c(^;w(t„,t„)) -$c(^;u(in)) - iv{tn,tn) - U{tr,)) 

^ .At ^ .At , ,,, ^ .At . 

*A(^;^B(^;w(tn,in))) -*c(^;^'(^n,i„)) 



(47) 



At 
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First we find that 



by using (41 1 (with v{tn,tn) as initial data). 
Introduce the function V — V{x,t) satisfying 



Then the very last line of (47) can be written as 

.At 



At 

tn + l/2 



(*c(^; ■)-i)ov{tn,tn)-{<i>c{Y' ■)-i)°i^itn) 

{Vt{a)-uti<7))da 



,1 + 1/2 ,1 



" + 1/2 ^1 

2 



{{V{a) + u{a)){V{a) - u{a)))^ - - u{a))xxx)da. 



Taking the ^ norm we find 
,At 



< 



< 



,1 + 1/2 ,1 



tn + l/2 



(- \\{V{a) + u{a)){V{a) - + H^^) - 

(if ( \\V{<j)\\h.^s + h(cT)||^.-s ) \\V{a) - + \\V{a) - u{a)\\ 



,, + 1/2 



\\V{a) - u{a)\\fjs-6 da. 



By the H stability of the KdV equation 

<K\\v{t„,tn)-u{t^)\\Hs-e < K^At, 

since by the arguments of Section 2.3 we have the estimate ||u'(t„)||^s_6 < ll^'^(^n)||/fs 
Kr.At. Therefore 



(*c(^; ■)-l)ov{tr,,tr,)-(^^c{^; ■)-l)ou{tn) 
Thus we have shown that 

(48) ||w(i„+i/2) - w(i,0L=-9 < Ka^At^. 

The other terms on the right hand side of (45 1 can be estimated in the same manner, 
using the expressions ( 46 ) . Thus we conclude that (cf . ( 44 1 ) 

(49) Wwrn^^^^KK^At". 

If t,r e [t„,t„+i/2], we have \\v{t,T)\\^k = \\v{t,t)\\^u, and if i,T e [t 

n+1/2 ; i 

then an estimate analogous to (16 1 shows that 

\Ht,T)\\H-,^\\v{t,t)\\H>\ <Ko.\t-T\. 

The rest of the argument follows the procedure for the Godunov splitting. Now 
s — 9 > fc is the same as s > 17. Assuming this, we get 



,At 



< Kc,Ar 



\\v{t,T)\\^-, < \\v{t,t)\\^, + \\vit,T)\\^, - \\vit,t)\ 



< K + A' At^ + Kr,At. 
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Choosing a such that K < a/4, and then At such that KaAt{At +1) < a/4 



implies that t)!]^:^;; < a/2. Hence by the bootstrap lemma and (44l, the 
following holds: 

Theorem 3.5. Fix T > 0. Let uq S for some s > 17. Then for At sufficiently 
small we have 

\\v{t,t)-u{t)\\j,,^o <KAt^, [0,T], 
where the constant K depends on Uq, s and T only. 
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